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It is shown that a first-order cosmological perturbation theory for the open, flat and closed 
Friedmann-Lemaitre- Robertson- Walker universes admits one, and only one, gauge-invariant variable 
which describes the perturbation to the energy density and which becomes equal to the usual 
Newtonian energy density in the non-relativistic limit. The same holds true for the perturbation 
to the particle number density. Using these two new variables, a new manifestly gauge-invariant 
cosmological perturbation theory has been developed. 

Density perturbations evolve diabatically. Perturbations in the total energy density are gravita- 
tionally coupled to perturbations in the particle number density, irrespective of the nature of the 
particles. There is, in first-order, no back-reaction of perturbations to the global expansion of the 
universe. 

Small-scale perturbations in the radiation-dominated era oscillate with an increasing amplitude, 
whereas in older, less precise treatments, oscillating perturbations are found with a decreasing 
amplitude. This is a completely new and, obviously, important result, since it makes it possible to 
explain and understand the formation of massive stars after decoupling of matter and radiation. 

PACS numbers: 04.25.Nx, 97.10.Bt, 98.62.Ai, 98.80. Jk 



I. INTRODUCTION 

Since measurements of the fundamental parameters of 
our universe are today very precise, cosmology is nowa- 
days a mature branch of astrophysics. Despite advances 
in observational as well as theoretical cosmology, there is, 
as yet, no manifestly covariant cosmological perturbation 
theory. It is the purpose of this article to fill up this gap. 

We redo the calculations of Lifshitz [1] and Lifshitz and 
Khalatnikov [2], and we use new results of the literature, 
which were not known at the time when Lifshitz and Kha- 
latnikov developed their theory. Combining these results 
with new insights obtained by us and described in detail 
in our report [3], we are able to develop a perturbation 
theory with the properties that both the evolution equa- 
tions and their solutions are invariant under general in- 
finitesimal coordinate transformations a;^ — {x'^) . 
We refer to such a theory as a manifestly covariant gauge- 
invariant perturbation theory. 

A. Former Insights 

Firstly, it is mandatory [4-9] to use gauge-invariant 
variables to construct a perturbation theory. Secondly, 
although Lifshitz and Khalatnikov were aware of the fact 
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that the system of perturbation equations can be divided 
into three independent systems of equations, namely for 
gravitational waves, vortices, and scalar perturbations, 
they did not use explicitly the decompositions (11) and 
(13) found by York [10] and Stewart [11]. This decom- 
position makes the computations much more tractable. 
Finally, an important result [10, 11] is that the perturbed 
metric tensor for scalar perturbations can be written in 
terms of two potentials (12). This facilitates the deriva- 
tion of the Newtonian results (56) and (57) in the non- 
relativistic limit. 



B. Our New Insights 

In addition to the results found in the literature, we 
use the new insights obtained by us [3] . Firstly, we have 
found, using the combined First and Second Law of Ther- 
modynamics, that, in general, density perturbations do 
not evolve adiabatically, as we will demonstrate in Sec- 
tion VII. Therefore, it is essential that in a cosmological 
perturbation theory one uses an equation of state for the 
pressure p = p{n, e) rather than p — p{e), where n is the 
particle number density and e is the energy density of the 
universe. This enables us to show that a small negative 
internal pressure fluctuation may have a favorable effect 
on the growth rate of a density perturbation. Secondly, we 
have found that for scalar perturbations, the perturbed 
momentum constraint equations (8b) can be rewritten 
as an evolution equation (14b) for the local perturbation 
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(9) to the global spatial curvature ^R(o) (6)- Using 
this new insight we can rewrite the evolution equations 
for scalar perturbations as simple extensions (14) of the 
background equations (5). The sets of equations (5) and 
(14) are written with respect to the same system of ref- 
erence. From these two sets of equations it follows that 
only three independent scalars, namely the energy den- 
sity e, the particle number density n, and the expansion 
scalar 0, play a role in a perturbation theory. This re- 
duces the number of possible gauge-invariant variables 
considerably. Finally, in a first-order perturbation the- 
ory the gravitational field is already weak. Therefore, we 
only need to take the limit that all particle velocities are 
negligible with respect to the speed of light, u/c ^ 0, the 
so-called non-relativistic limit, see (51), in order to show 
that our treatment reduces to the Newtonian theory of 
gravity. In fact, we show that there exist two, and only 
two, gauge-invariant variables defined by (17a) with the 
property that in the non-relativistic limit the well-known 
Newtonian results (56) and (57) show up. Consequently, 
efi') and nfl^ are the real, physical, energy density pertur- 
bation and the particle number density perturbation. 



II. RESULTS 

The results of our approach are as follows. Pertur- 
bations in the total energy density are gravitationally 
coupled to perturbations in the particle number density 
(21b), irrespective of the nature of the particles [i.e., or- 
dinary matter or Cold Dark Matter (cdm)] and inde- 
pendent of the scale of the perturbations. Local density 
perturbations do not affect the global expansion of the 
universe: in first-order there is no back-reaction. 

For large-scale perturbations, our new theory corrob- 
orates the outcomes of the standard theory, with the ex- 
ception of the fact that we do not find the non-physical 
gauge-modes: see our comments below expressions (34) 
and (47). 

For small-scale perturbations, however, the matter is 
quite different. In the literature it is assumed that in the 
small-scale limit the Newtonian theory of gravity applies 
and that gauge dependent variables become in the small- 
scale limit automatically gauge-independent [4] . This as- 
sumption is not true. Indeed, we show in Section IX 
that a gauge dependent quantity is still gauge depen- 
dent in the non-relativistic limit. Consequently, a mani- 
festly gauge-invariant theory is of the utmost importance 
for the study of the evolution of small-scale perturba- 
tions. For small-scale perturbations our perturbation 
theory yields results that differ substantially from the 
standard theory. Firstly, in the radiation-dominated era 
of a flat FLRW universe we find that small-scale perturba- 
tions oscillate with an increasing amplitude according to 
(35). The outcomes of the standard theory, namely de- 
caying perturbations, are found to be a coordinate arte- 
fact. Secondly, after decoupling of matter and radiation, 
small-scale density perturbations evolve also very differ- 



ently from the behavior predicted by the standard theory. 
For small-scale perturbations, our new evolution equa- 
tion (37a) contains an entropy term which is of the same 
order of magnitude as the pressure term. This will shed 
new light on the problem of structure formation in the 
universe. 



III. BASIC EQUATIONS IN SYNCHRONOUS 
COORDINATES 

Due to the general covariance of the Einstein equa- 
tions and conservation laws, Einstein's gravitational the- 
ory is invariant under a general coordinate transforma- 
tion a;'* — >■ x^{x'^), implying that preferred coordinate 
systems do not exist. In particular, the linearized Ein- 
stein equations and conservation laws are invariant under 
a general linear coordinate transformation 

x° x° -^p{t,x), x' x' -C{t,x), (1) 

where 'ip{t,x) and ^'(t, a;) are four arbitrary, first-order 
('infinitesimal') functions of the time and space coordi- 
nates. Since preferred systems of reference do not exist 
and since our result (21) is manifestly gauge-invariant, 
one may use any suitable and convenient coordinate sys- 
tem to perform the calculations. In the Newtonian the- 
ory of gravity all systems of reference are synchronous 
because space and time are decoupled in this theory, see 
(54). In order to facilitate the interpretation of our new 
gauge-invariant quantities, we use, just like Lifshitz and 
Khalatnikov do, synchronous coordinates. In this coor- 
dinate system the metric of flrw universes has the form 

300 = 1, 5oj = 0, Qij ^ -a^{t)g,j{x), (2) 

where a{t) is the scale factor of the universe, and gij is 
the metric of the three-dimensional maximally symmetric 
subspaces of constant time. The functions a-^d C of 
the general infinitesimal transformation (1) reduce to 

^P = ^P[x), e = g^'dki^ix) ^^7^ + (3) 

J a^iV 

if only transformations between synchronous coordinates 
are allowed. In (3), ip{x) and x'(a;) are four arbitrary, 
infinitesimal, functions of the spatial coordinates. 



A. Background Equations 

The complete set of zeroth-order Einstein equations 
and conservation laws for an open, flat or closed 
FLRW universe filled with a perfect fluid with energy- 
momentum tensor 

Ti''' = {e + p)u''u'' - pg""" , p = p{n,e), (4) 
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is, in synchronous coordinates, given by 



Constraint: 
Evolution: 
Conservation: 
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2 -"-(0) 



K6 n 



■A, 



£(0) = -3i7e(o)(l + w), 
72(0) — "3i?n(o) , 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 



where the equation of state for the pressure p is sup- 
posed to be a given function of the energy density e 
and the particle number density n. The Goi constraint 
equations and the Gij, i ^ j, dynamical equations are 
identically satisfied. The Ga dynamical equations are 
equivalent to the time-derivative of (5a). Therefore, 
the Gij dynamical equations need not be taken into ac- 
count. In equations (5) A is the cosmological constant. 



K = 8ttG/c*, and 



P(o)/£(o)- An overdot denotes 



differentiation with respect to ct and the sub-index (0) 
refers to the background, i.e., unperturbed, quantities. 
Furthermore, H = a/a is the Hubble function which is 
equal to H = |0(o), where 9^o) is the background value of 
the expansion scalar 6* = u^.^ with the four- velocity, 
normalized to unity (u^u^ = 1)- A semicolon denotes 
covariant differentiation with respect to the background 
metric (?(o)/ii/. The spatial curvature ^i?(o) is given by 



- 

^(0) — 



6fc 



k = -1,0, +1. 



(6) 



The quantity ??(o) is the three-divergence of the spatial 
part of the four- vector u'^g^ . For an isotropically expand- 
ing universe we have m^o) = 6^o, so that d^o) = 0. 

From the system (5) one may infer that the evolution of 
an unperturbed flrw universe is determined by exactly 
three independent scalars, namely 



'Uf,u^, n^N^'u^, e^u''.^, (7) 

where = nu^^ is the cosmological particle current 
four-vector, which satisfies the particle number conser- 



vation law N^- 



0. 



B. Perturbation Equations 

The complete set of first-order Einstein equations and 
conservation laws for the open, flat and closed flrw uni- 
verse is, in synchronous coordinates, given by 



k\i ~ i\k — 2k(£(o) +P(0))W(i)i, 

£(1) + 3i?(e(i) + Pd)) + (£(0) + P(o))Oa) = 0, 
1 d r ^1 
c di °' +P(o))"(i) ~ 

9lo)P(i)\k + 5i?(£(o) +Pfo))u\,) = 0, 

71(1, + 3ifn(i) + n(o)6'(i) = 0, 



(8a 
(8b 

(8c 
(8d 



(8f: 



where h'^j = glo^hkj is the perturbed metric, and g^^^ 
the background metric (2) for an open, closed or flat 
FLRW universe. Quantities with a sub-index (1) are the 
first-order counterparts of the corresponding background 
quantities with a sub-index (0). The first-order pertur- 
bation to the spatial part of the Ricci tensor is ^Rl^j and 
its trace is given by 



— 9li){h^k\i\j — h^i\j\k) + I ^R{0)h^k 



(9) 



A vertical bar denotes covariant differentiation with re- 
spect to the background metric g^onj- Expression (9) 
is the local perturbation to the global spatial curvature 
'^R(o) due to a local density perturbation. Finally, is 
the first-order perturbation to the expansion scalar 9. It 
is found that 



1 Uk^ 



*(l)|fc' 



(10) 



where 79(i) is the divergence of the spatial part of the 
perturbed four- vector u^,. The quantities (9) and (10) 
play an important role in the derivation of our manifestly 
gauge-invariant perturbation theory. 



IV. FIRST-ORDER EQUATIONS FOR SCALAR 
PERTURBATIONS 



Since the work of Lifshitz and Khalatnikov, it is known 
that for FLRW universes, the set (8) can be broken down 
into three independent systems of equations, namely a 
system for gravitational waves (tensor perturbations), for 
vortices (vector perturbations), and for scalar perturba- 
tions. In addition they showed that only scalar pertur- 
bations are coupled to density perturbations. Lifshitz 
and Khalatnikov used spherical harmonics to classify the 
three types of perturbations. A disadvantage of their 
approach is that using spherical harmonics is computa- 
tionally demanding. A new insight, gained in 1974 by 
York [10], and in 1990 elucidated by Stewart [11], makes 
the computations much more tractable. York and Stew- 
art showed that the symmetric perturbation tensor h^j 
can uniquely be decomposed into three parts, i.e., 



^±k = ^*fe = 



0, 



h^k = 0, 



(11a) 
(lib) 



where the scalar, vectorial and tensorial perturbations 
are denoted by ||, _L and *, respectively. Moreover, they 
demonstrated that the component h'^y can be written in 
terms of two independent potentials 0(t, x) and C(i, a;): 



hi 



(12) 



Finally, Stewart also proved that the spatial part of the 
perturbed four-vector uf-^^ can uniquely be decomposed 
into two parts 

M(i) = M(i)|| + M(i)_L, (13a) 
V • = V • , V X M(i, = V X ■U(i)_L, (13b) 
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where (V/)* = g'^djf. 

The perturbed Ricci tensor, ^R\^j, being a symmetric 
tensor, should also obey the decomposition (11a) with 
the properties (lib), i.e., ^R^^^j^k = As a consequence, 
/ij^ must obey /i^'ifcii ~ addition to (lib), as follows 
from (9). This additional property is a necessary and suf- 
ficient condition for the decomposition of the perturba- 
tion equations into three distinct types of perturbations, 
with coupled to scalar perturbations and M(i)J_ cou- 
pled to vortices, according to their properties (13b). 

Since only the scalar part of the perturbations is cou- 
pled to density perturbations, we may replace in (8)-(10) 
by /ly^ and m'^j by to obtain perturbation equa- 
tions which exclusively describe the evolution of scalar 
perturbations. Using the decompositions (11a) and (13a) 
and the properties (lib) and (13b), we can rewrite the 
evolution equations for scalar perturbations in the form 



(14a) 



2-H^(^'(i) - - \ ^R{i)\\ — i^>^(i), 
^R,,,\\+2H%,,\\- 

2k£(o)(1 + w)^^,, + I %oAdii) - ^(1)) = 0, (14b) 

£(1, + 3i?(£(i) + P(i,) + £(0,(1 + W)e^r> = 0, (14c) 

^<,+i/(2-3^^).,, + ^-^^=0, (14d) 

+ 3iJn(i) n(o)6'(i) = 0. (14e) 

The set (14), which is the perturbed counterpart of 
the set (5), consists of five equations for the five un- 
known functions £(i,, ri(i), and 9^ly The first- 
order perturbation to the pressure is given by the per- 
turbed equation of state = Pnn^i) + Pe£(i), where 
Pn = {dp/dn)f; and p^ = {dp/de)n are the partial deriva- 
tives of the equation of state p{n, £), and = p^o)l£(o)- 
Finally, the symbol denotes the generalized Laplace 
operator with respect to the three-space metric gij. 

We now sketch the derivation of the basic equations 
(14) for scalar perturbations. Eliminating h^^f, from (8a) 
with the help of (10) yields (14a). Multiplying both sides 
of equation (8b) by g\l^ and taking the covariant deriva- 
tive with respect to the index j and using (9), (10) and 
the fact that the background connection coefficients Y^^ij 
are for flrw metrics independent of time, we arrive at 
equation (14b). In the derivation of (14b) we have used 
that g\i)h^i\j\k = g^h'' i\k\j , which follows from g'^l^ |j. = 
and the symmetry of K^^ . Equations (8c) need not be 
considered, since for i ^ j they are not coupled to scalar 
perturbations, whereas the trace of (8c) is, just as in the 
background case, equivalent to the time-derivative of the 
constraint equation (14a). Finally, equation (14d) can 
be obtained from (8c) by taking the covariant deriva- 
tive with respect to the metric g{a)ij- This concludes the 
derivation of the system (14). As follows from its deriva- 
tion, this system is, for scalar perturbations, equivalent 
to the full set of first-order Einstein equations and con- 
servation laws (8). 



UNIQUE COSMOLOGICAL DENSITY 
PERTURBATIONS 



The background equations (5) and the perturbation 
equations (14) are both written with respect to the same 
system of reference. Therefore, these two sets can be 
combined to describe the evolution of the five background 
quantities 6*10) = 3i/, ^R^q), £(o), 'd{o) — and n^g^, and 
their first-order counterparts 9^^), '^(i) and 

ri(i). Just as in the background case, we again come 
across the three independent scalars (7). Consequently, 
the evolution of cosmological density perturbations is de- 
scribed by the three independent scalars (7). A compli- 
cating factor is that the first-order quantities £(i) and 
71(1), which are supposed to describe the energy density 
and the particle number density perturbations, have no 
physical significance, as we will now establish. 

A first-order perturbation to one of the scalars (7) 
transforms under a general (not necessarily between syn- 
chronous coordinates) infinitesimal coordinate transfor- 
mation (1) as 



(15) 



where S^o) and S^l-| are the background and first-order 
perturbation of one of the three scalars S = e,n,9. In 
(15) the term S = tpS^o) is the so-called gauge mode. The 
complete set of gauge modes is given by 



#(0), 



(1) — V"^(0)7 

'Ad II =4-ff 



'(1) — V'^CO); 



(16a) 
(16b) 



where expressions (16b) hold true only in synchronous 
coordinates. The quantities (16) are mere coordinate 
artifacts, which have no physical meaning, since the 
gauge function 'ijj{x) is an arbitrary (infinitesimal) func- 
tion. Equations (14) are invariant under coordinate 
transformations (1) with (3). This property combined 
with the linearity of the perturbation equations, im- 
plies that a solution set (£(i), n(i), 0(i), i?(i), ^i?(i)||) can 
be augmented with the corresponding gauge modes (16) 
to obtain a new solution set. Therefore, the solution 
set (£(1), 0(1), ??(i), '^i?(i)|| ) has no physical significance, 
since the general solution of the set (14) can be modified 
by an infinitesimal coordinate transformation. 

It is possible to eliminate the gauge modes for the 
scalars £, n and 9, by constructing so-called gauge- 
invariant variables, i.e., variables that do not change un- 
der general infinitesimal coordinate transformations (1). 
All attempts [4-9] to construct a gauge-invariant cos- 
mological perturbation theory yield different outcomes. 
For instance, Bardeen [4] defines three different gauge- 
invariant variables to describe energy density perturba- 
tions, which become equal to each other only in the small- 
scale limit. The variables of Bardeen differ from the ones 
used by Mukhanov et al. [5, 9], and the variables used 
by Mukhanov et al. differ, in turn, from the ones used by 
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Ellis et al. [6-8] . Apparently, the number of possibilities 
to construct gauge-invariant variables is large. This is 
why there is no consensus about which variable exactly 
describes the evolution of energy density perturbations 
in the universe. This is the notorious gauge problem of 
cosmology. 

To unravel this problem, we first have rewritten the 
perturbation equations (8) into the form (14) in or- 
der to isolate the scalar perturbations from the vor- 
tices and gravitational waves. This reduces the num- 
ber of possible gauge-invariant variables substantially, 
since we need only consider the three independent scalars 
(7). All scalar perturbations transform in the same way 
(15). This implies that we can combine two independent 
scalars to eliminate the gauge function ip{t, x). With the 
three independent scalars (7), we can make three differ- 
ent sets of three gauge-invariant variables. In each of 
these sets exactly one gauge-invariant quantity vanishes. 
As we will show in Section IX, the only set for which 
the corresponding perturbation theory yields in the non- 
relativistic limit u/c — >■ the well-known Newtonian re- 
sults (56) and (57) is given by 



fc(i) 



£(1) 



7(0) 



^(1) 



0. 



(17a) 



(17b) 



7(0) 



It follows from the general transformation rule (15) that 
the quantities (17) are invariant under the infinitesimal 
transformation (1), i.e., they are gauge-invariant, hence 
the superscript 'gi'. The physical interpretation of (17b) 
is that, in first-order, the global expansion O^oy — 3H is 
not affected by a local density perturbation. 

The gauge-invariant quantities (17a) are completely 
determined by the background equations (5) and their 
perturbed counterparts (14). In principle, we can use 
these two sets to study the evolution of density perturba- 
tions in FLRW universes. The set (14) is still too compli- 
cated, since it also admits the non- physical solutions (16). 
Our aim will be a system of equations for ef') and nfl^ 
that do not have the gauge modes (16) as solution. In 
other words, we are looking for a perturbation theory for 
which not only the differential equations are invariant un- 
der general infinitesimal coordinate transformations (1), 
but also their solutions. We refer to such a theory as 
a manifestly covariant gauge-invariant perturbation the- 
ory. The construction of such a theory will be the subject 
of the next section. 



VI. PERTURBATION EQUATIONS FOR 
DENSITY PERTURBATIONS 

In this section we will outline the derivation of our new 
perturbation theory. Firstly, we observe that the gauge 
dependent variable 0(i) is not needed in our calculations, 
since its gauge-invariant counterpart 9fl^, (l^b), vanishes 



identically. Eliminating from the differential equa- 
tions (14b)-(14e) with the help of the algebraic equation 
(14a) yields the set of four first-order differential equa- 
tions 



2Ke(o)(l + w)i9(i) + 



-n-(o) 



3H 



ine. 



Pi?<,),l) =0, (18a) 



£(1) +3i7(£(i) +P(i)) + 



£(0,(1 + w) l?(i) 



1 

2H 



(k£(i, 



+ i/(2 - 3/32)^(,, -I- 



1 3 

2 -''-(1) 
72 



'R(i 
V'P(i) 



-0, 



)(l + u;) 



a' 



'^(i) + ^ ('^^(1) + 5 ^^(i)ll) 



0, 



(18b) 
(18c) 

(18d) 



for the four quantities £(i,, n^,, z^^j and ^-R(i)||. 

Using the background equations (5) to eliminate all 
time-derivatives and the first-order constraint equation 
(14a) to eliminate we can rewrite the gauge- invariant 
quantities (17a) as 



- 3£(o)(l -f w){2Hd,,, + \ 



+3k£(o)(1 + w) 

371(Q) ('t£(i) 



^i?(0) -I- 3k£(o)(1 + w) 



(19a) 
(19b) 



These quantities are now completely determined by 
the background equations (5) and the first-order equa- 
tions (18). In the study of the evolution of density pertur- 
bations, it is convenient not to use efl^ and nfl^ directly, 
but instead their corresponding contrast functions 6^ and 
dn defined by 



d^{t,x) 



£(1) (^) ''') 



"(O)(0 



(20) 



We now rewrite the system of equations (18) for the four 
independent quantities £(i,, n^iy, and ■^i?(i)|| into a new 
system of equations for the two independent quantities 
and S„. The final result is our new manifestly covariant 
gauge-invariant perturbation theory: 



biK + b25e = 63 



Sn - 



1 



1 d 
cdt 



Sn - 



1 



£(0)(1 + w) 



Sn - 



1 



(21a) 
(21b) 



These are two differential equations for the two indepen- 
dent and gauge-invariant quantities 5^ and Sn- It follows 
from equation (21b) that perturbations in the total en- 
ergy density are gravitationally coupled to perturbations 
in the particle number density if p„ = (dp/dn)^ < 0. 
This is the case in a flrw universe in the radiation- 
dominated era and after decoupling of matter and radi- 
ation. This coupling is independent of the nature of the 
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particles, i.e., it holds true for ordinary matter as well as 

CDM. 

The system of equations (21) is equivalent to a sys- 
tem of three first-order differential equations, whereas the 
original set (18) is a fourth-OTder system. This difference 



is due to the fact that the gauge modes, which are so- 
lutions of the set (18), are completely removed from the 
solution set of (21): one degree of freedom, namely the 
gauge function ip, has disappeared. 

The coefficients &i, 62 and 63 occurring in equation 
(21a) are given by 



1 



2H{1 + 3/32) 



^3H 3^(0) +3Ke(o)(l + w) 
&2 = -5Ke(o)(l + w){l + 3w) + {l~3w + 6p^{2 + 3w)) 



■6H^ [ w + 



K£(o)(l + W) 



b. = 



^i?(0) + 3/t£(o)(l + w) 



%o, ( '^w + 



i72(i + 6w)(l + 3/32) 
^i?(o) -t- 3Ke(o)(l + w) 



37? 



3Ke(o)(l + w) 



2pn P 



(22a) 



(22b) 





"(0) / 


+ Pn 






"£(0) I 



V2 

d 



2 -pi?(o, , (22c) 



where p„„ = d'^p/dv? and p^n = d'^p/de dn. 

The background equations (5) and the new perturba- 
tion equations (21) constitute a set of equations which 
enables us to study the evolution of small fluctuations in 
an open, flat or closed flrw universe filled with a perfect 
fluid with an equation of state p(n, e). 



VII. DIABATIC DENSITY PERTURBATIONS 

In this section we will relate equations (21) to ther- 
modynamics and we will show that, in general, density 
perturbations do not evolve adiabatically, as is usually 
assumed, but diabatically, i.e., they exchange heat with 
their environment during their evolution. 

The combined first and second law of thermodynamics 
is given by 



TAs 



(23) 



where s is the entropy per particle. Using the background 
equations (5) it follows that S(o) — 0. This implies with 
(15), that S(i) = s^i) is automatically gauge-invariant. 
The thermodynamic relation (23) can, using (20), be 
rewritten in the form 



1 



(24) 



Thus, the right-hand side of (21a) is related to local per- 
turbations in the entropy. 

Adiabatic perturbations do not exchange heat with 
their surroundings, so that T^o-^sf^^ — 0. Using (20) and 
the background conservation laws (5c) and (5e), we find 



from 



the adiabatic condition 



'(0)^(1) 



0. 



(25) 



In a non-static universe (i.e., £(0, ^ and rijo) ^ 0) this 
equation is fulfilled if, and only if, the energy density is a 
function of the particle number density only, i.e., if and 
only if e = e(n). This can be demonstrated as follows. 
From thermodynamics it is known that e = e{n, T) and 
p = p{n, T), where the particle number density n and the 
temperature T are independent quantities. Substituting 
e — e{n,T) into (25) yields 



"'(o)-'(i) ''(i)-f(o 



= 0. 



(26) 



Since n and T are independent quantities, the adiabatic 
condition (26) is satisfied if, and only if. 



-] =0 



(27) 



In particular, in the non- 
and p = 0, density per- 



implying that s = e{n). 
relativistic limit, where e — nnic 
turbations are adiabatic. In all other cases, e — e{n, T) 
and p — p{n, T) [hence p — p{n, e)], local density pertur- 
bations evolve diabatically. 



VIII. THE FLAT FLRW UNIVERSE 

We consider a fiat [k — 0, (6)] flrw universe in its 
radiation-dominated phase and in the era after decou- 
pling of matter and radiation. In this section only, we 
take A = 0. 



Radiation-dominated Era 



In this epoch we have e 
body constant and the radiation temperature. 



O'bT^j where ae is the black 



The 
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pressure is p = ^e, so that p„ = 0, Pe = ^ and = |. 
In this case, the perturbation equations (21) reduce to 



1 d_ 
cdt 



) (I 



IS. 



s, = o, 

0. 



(28a) 
(28b) 



Equation (28b) implies that the difference (5„ — '^S^ de- 
pends only on the spatial coordinates. Since the universe 
is non-static and {de/dT)n 0, it follows that (26) can- 
not be satisfied, implying that the perturbations are not 
adiabatic, so that Sn{to,x) — ^6e{to,x) ^ 0. 

Substituting (5e (i, x) = Sg{t,q)exp{iq-x) into equation 
(28a) and using the well-known solutions of the back- 
ground equations (5) 



H (xf- 



oc t 



(X t 



-3/2 



equation (28a) reduces to 



where /ir is given by 

27r 1 

Mr 



At 



1 



1 

2^2 



a (X t 



0, 



1/2 



Xq = Xa{to), 



(29) 



(30) 



(31) 



with Xa(to) the physical scale of a perturbation at time 
to, and |q| = 2tt/X. The exact solution of (30) is 

Se{T, q) = Ai{q) sin {fi^^/r) + ^2(9) cos (/ir%/r) ^/T, 

(32) 

where r = t/to- The 'constants' of integration Ai{q) and 
j42(q) are given by 



^1,2(9) = 5,{to,q) 



sm /Zr 
cos /ir 



1 cos /if 
//r Sin/ir 



Seih.q) 



Se{tQ,q) 

H{to) 



(33) 



For large-scale perturbations (A — 00), we arrive at 

Seito) 



2Seito) 



H{to) 

Sejto) 
H{to) 



t 
to 



(34) 



The energy density contrast has two contributions to the 
growth rate, one proportional to t and one proportional 
to These two solutions have been found, with the 

exception of the precise factors of proportionality, by a 
large number of authors. See Lifshitz and Khalatnikov 
[2], (8.11), Adams and Canute [12], (4.5b), Olson [13], 
page 329, Peebles [14], (86.20), Kolb and Turner [15], 
(9.121) and Press and Vishniac [16], (33). 



A new result of our perturbation theory is that small- 
scale perturbations oscillate with an increasing ampli- 
tude, proportional to i^/^, according to 



Se{t,q) W (Se(to,q) 



Mr - Mr 



(35) 



These small-scale perturbations manifest themselves as 
small temperature fluctuations in the cosmic background 
radiation. 

In contrast to our theory, the standard perturbation 
theory predicts oscillating density perturbations with a 
decreasing amplitude. As is explained in our report [3], a 
decreasing amplitude is due to the spurious gauge modes 
present in the standard theory. This is why a manifestly 
gauge-invariant perturbation theory is of the utmost im- 
portance. 



B. Era after Decoupling 

In this section it is assumed that the CDM particle mass 
is larger than or equal to the proton mass, tocdm > 'tih, 
implying that for the mean particle mass m we have 
mc^ > k-QT. 

Once protons and electrons combine to yield hydro- 
gen, the radiation pressure becomes negligible, and the 
equations of state become those of a non-relativistic 
monatomic perfect gas 

e{n,T) = nmc^ + ^nkBT, p{n,T) ~ nk^T, (36) 

where fee is Boltzmann's constant, m the mean particle 
mass, and T the temperature of the matter. We have 
= I and pn — ~^mc'^, so that the system (21) re- 
duces to 



2 V2 

3 a 



2 (Sn - Se), 



1 d_ 
cdt 



(Sn - S,) = -2H (Sn - S,) , 



(37a) 
(37b) 



where (3 = \/P{a)/^{o) is, to a good approximation, given 

by 



m 



Vs{t) 



5 kBTfo^jt) 

3 TOC^ ' 



T,o, cx a 



(38) 



with I's the adiabatic speed of sound and r(o) the matter 
temperature. Equation (37b) implies 



Sn - Se 



(39) 



where we have used that H = a /a. From the expressions 
(36) it follows that the relative pressure perturbation is 
given by 



5n + 5t, 



(40) 
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where 5t is the relative matter temperature perturbation 
defined by 5t = TfljT^o) and Tfl^ is defined as ef', in 
(17a) with e replaced by T. Using that after decoupling 
^bT(o)/(»tic^) ^ 1, we find from (36) the relation 



5„ 



3 feBTfo) 



5t- 



(41) 



2 mc^ 

Combining (38) and (39) we find from (41) that 

5T{t,x)^5T{h,x), (42) 

to a very good approximation. Using the well-known 
solutions of the background equations (5) 



H cat ^ , e,o, cx t 



71(0) OC t 



a oc t 



2/3 



(43) 



and 6{t,x) = S{t,q) exp{iq ■ x), equations (37) can be 
combined into one equation 



T 



4 



9^8/3 



10 

9^2 



<5e = -^^^T(^o,g), (44) 



where t = t/tQ and a prime denotes differentiation with 
respect to t. The parameter is given by 



IJ-n 



2tt 1 ws(to) 



An 



Ao H{to) c 

The general solution of equation (44) is found to be 
5,{T,q) = [Bi(q)J 7(2M,nr-i/3) 



(45) 



2 

5^2/3 



-1/3^ 



^1/2 



6Tito,q)i 



(46) 



where ^±7/2(2;) are Bessel functions of the first kind and 
Bi{q) and B2{q) are the 'constants' of integration. 

In the large-scale limit A — >■ 00 terms with vanish, 
so that the general solution of equation (44) is 



H{to) 



2 

t \ 3 



(47) 



Thus, for large-scale perturbations the initial value 
'^t(^Oj Q) does not play a role during the evolution: large- 
scale perturbations evolve only under the influence of 
gravity: they are so large that neither heat exchange 
nor pressure perturbations do play a role during their 
evolution in the linear phase. For perturbations much 
larger than the Jeans scale, gravity alone is insufficient 
to account for structure formation within 13.75 Gyr. The 
solution proportional to t^^^ is a standard result. Since 
(5e is gauge-invariant, the standard non-physical gauge 



mode proportional to is absent from our theory. In- 
stead, a physical mode proportional to t~^^^ is found. 
This mode has also been found by Bardeen [4], Table I, 
and by Mukhanov et al. [5], expression (5.33). In or- 
der to arrive at the t~^/^ mode, Bardeen has to use the 
'uniform Hubble constant gauge.' In our treatment the 
Hubble function is automatically uniform, without any 
additional gauge condition, see (17b). 
In the small-scale limit A ^ 0, we find 



5e{t,q)~-UT{to,q) 



1 

'3 r 



l5T{h,q)+5,{to,q) 



n 



(48) 



Thus, density perturbations with scales much smaller 
than the Jeans scale oscillate with a decaying amplitude 
which is smaller than unity: they are so small that grav- 
ity is insufficient to let perturbations grow. Pressure per- 
turbations alone cannot make these perturbations grow. 
Consequently, small-scale perturbations will never reach 
the non-linear regime. 

The study of the evolution of perturbations with in- 
termediate scales will be conducted in a forthcoming pa- 
per by solving equation (44) numerically. It is demon- 
strated that for density perturbations with scales some- 
what larger than the Jeans scale, the action of both grav- 
ity and an initially small negative pressure perturbation 
together may result in massive stars several hundred mil- 
lion years after decoupling of matter and radiation. 



IX. NON-RELATIVISTIC LIMIT 

In Section V we have shown that only two gauge- 
invariant quantities efl-^ and v!fl) exist, which could be 
the real energy density and particle number density per- 
turbations. In this section we will demonstrate that in 
the non-relativistic limit u/c — >■ our theory reduces to 
the well-known results (56) and (57) of the Newtonian 
theory of gravity and that the quantities efl-^ and nf', 
become equal to their Newtonian counterparts. This im- 
plies that efi) and nf^^ are indeed the local perturbations 
to the energy density and particle number density per- 
turbations. 

We use the result (12) found by Stewart [11], which 
implies that, in general, cosmological density perturba- 
tions are described by two potentials </> and C,. In the 
non-relativistic limit of a flat {k — 0) flrw universe, the 
potential C disappears from our theory, thus paving the 
way towards the Newtonian theory of gravity. 

As follows from (9) and (10), both potentials occur 
only in the quantities ^R{i)\\ and O^^y A first step to get 
rid of the potential C is to consider a flat (k = 0) flrw 
universe. For a fiat universe, we find that, using (12), the 



9 



local perturbation to the spatial curvature (9) reduces to 

4 V^(f> 



^7? 



(49) 



where is the usual Laplace operator. Using this ex- 
pression, the perturbation equations (14) read, for a flat 
FLRW universe. 



£(1) + 3i?(e(i 



AttG 



1 



47rG 



2 £(0,(1 + w)l?(i) = 0, 

P(i)) +e(o)(l + 'i«)6'(i) 

2 



+ i/(2 - 3/3^)z9(i, 
«(i) + 3iJn(i, + n(o)6'( 



0, 



1 V2p,„ 



£(o)(l + u;) a2 
= 0, 



0, 



(50a) 

(50b) 
(50c) 

(50d) 

(50e) 



where we have used (17a) to eliminate Sf^ from the con- 
straint equation (14a). In these equations, the poten- 
tial C occurs only in the quantity 6^^^, see (10) and (12). 
This potential disappears from our equations in the non- 
relativistic limit, as we will now show. The spatial part 
m^j^jII of the fluid four- velocity is gauge dependent with a 
physical component and a non-physical gauge part. We 
define the non-relativistic limit w/c — >■ by 



(1)11 physical 



r/' „ 

(1) II physical 



0, 



(51) 



i.e., the physical part of the spatial fluid velocity is neg- 
ligible with respect to the speed of light. In this limit, 
the kinetic energy per particle ^m{v'^) — '^k^T — > is 
very small compared to the rest energy mc^ per parti- 
cle, implying that the pressure p = nk-oT — > (n 7^ 0) 
is also vanishingly small. Substituting p — Q into the 
momentum conservation law (8e) yields, using also the 
background equation (5c) with w = P(o)/£(o) — 0, 



-2Hu 



(Dil- 



(52) 



Since the physical part of li^^jn vanishes in the non- 
relativistic limit, the general solution of equation (52) 
is exactly equal to the gauge mode 



(i)llgaugc 



(53) 



where we have used that H = d/a. Thus, in the non- 
relativistic limit (51) we are left with the gauge mode (53) 
only. Consequently, we may, without losing any physical 
information, put u^^_^ 



) II gauge 



equal to zero, implying that 
dktp = 0, so that is a constant in the non-relativistic 
limit. Therefore, the relativistic transformation (1) with 
(3) between synchronous coordinates reduces in the non- 
relativistic limit to the most general transformation 



(54) 



which is possible in the Newtonian theory of gravity. In 
(54), ip is an arbitrary constant and x*(^) three ar- 
bitrary functions of the spatial coordinates. 

Substituting i9(i) = (ufi,||)|fe = and p = into the 
system (50), we get in the non-relativistic limit 



Constraint: V-g 



AttG 



-a £?i) 



Evolution: V^(/) = 0, 
Conservation: £(i) + 3i/£(i) -I- £(o)fi'(i) = 0, 
71(1) + 3i7n(i) + n^o)0^r) = 0. 



(55a) 

(55b) 
(55c) 
(55d) 



The constraint equation (55a) can be found by subtract- 
ing times the time-derivative of the background 
constraint equation (5a) with ^R^o) = from the con- 
straint equation (50a) and using that 0(o) = 37? . 

Defining the potential ip{x) = (t){x) / {to) , equations 
(55a) and (55b) can be combined to yield 



efl{to,x) 



(56) 



which is the well-known Poisson equation of the 
Newtonian theory of gravity. 

Equations (55c) and (55d) have no physical significance 
since £(!,, n(i) and d^i^ are gauge dependent also in the 
non-relativistic limit, as follows from (5), (16a) and (54). 
Moreover, the physical equations (55a) and (55b) are de- 
coupled from the non-physical equations (55c) and (55d). 
Therefore, the latter two equations are not part of the 
Newtonian theory of gravity and need not be considered. 
Thus, in the non-relativistic limit, the potential C occur- 
ring in drops from our perturbation theory and we 
are left with one potential ip{x) = (j>{x) / (Iq) only. 

The expression (19a) reduces in the non-relativistic 
limit to £^) — — /(2k), which is, with (49), equiv- 
alent to the Poisson equation (56). Expression (19b) re- 
duces in the non-relativistic limit to the familiar result 



(57) 



where we have used that £(i) = n^y^mc^ and £(, 



n,o-.nic 



Finally, (55a) and (55b) imply that efl^ oc a~^. From 
(5c) with ui = we have £(o) cx a^^, so that with (20) 
and (57) we arrive at the well-known result (5^ = (5„ cx a. 

We have shown that our theory based on the gauge- 
invariant quantities (17) reduces in the non-relativistic 
limit d/c — to the well-known Newtonian results (56) 
and (57). Consequently, efl^ and nfl-, are the real, phys- 
ical, perturbations to the energy density and particle 
number density, respectively. 



X. STANDARD THEORY 

In order to compare the standard theory with the the- 
ory developed in this article, we consider a flat [k = 0, see 



(6)] FLRW universe after decoupling of matter and radi- 
ation. Since fcerjo) ^ 'tic^ the standard theory does not 
distinguish between perturbations in the particle number 
density and the total energy density. Therefore, we take 
in our perturbation theory (5„ — S^, implying that = 
and Sp ~ Sn, see (40) and (41). Equation (37a) reads 



0, 



(58) 



whereas (37b) is identically satisfied. This equation 
should be compared with the standard equation 



,5 + 2H6 



0. 



(59) 



Apart from the immediately obvious differences, namely 
the coefficients of H and ke^q-^ and the very small but 
unimportant difference between /3 and Ws/c, there are two 
major differences between (58) and (59). Firstly, (58) is 
derived from the full set of perturbed Einstein equations 
and conservation laws, without any assumptions or ap- 
proximations. In contrast to (58), the standard equation 
(59) can not be derived from the General Theory of Rel- 
ativity; instead, it is derived from the Newtonian theory 
of gravity, using dubious mathematics [17]. Secondly, the 
general solution of our equation (58) is gauge-invariant, 
whereas the general solution of the standard equation 
(59) depends on the gauge. This can be seen as follows. 
For large-scale perturbations, \7'^5 -H- 0, (59) reduces to 



(60) 



In contrast to (59), equation (60) can be derived from 
Einstein's gravitation theory: substituting 5 = e^i je^a^^ 
^(1) = 0, 71(0) = 71(1) = 0, P(o) = = 0, w ^ and 
'^R(o) = into the set (18) and using the background 
equations (5) yields (60). As a consequence, (60) is a 
relativistic equation. Using (43), we find for its general 
solution 



6{t) = Ct^/^ - 3H{to)^pT- 



(61) 



where C is an arbitrary constant, and r = t/tQ. Since 
the solutions of the set (18) are gauge dependent, the 
solution (61) of equation (60) is also gauge dependent 
and the mode proportional to is the gauge mode as 
follows from (16a) and £(0) oc t~'^. Note that — 
implies with (16b) that the gauge function ip is constant. 
Since the solution (61) can be obtained from the general 
solution of (59) by taking the large-scale limit, it follows 
that also the general solution of (59) contains the gauge 
function ip. Consequently, one cannot impose initial con- 
ditions, since the gauge function ijj cannot be fixed. We 
have to conclude that also the standard equation (59) 
has no physical significance and should, therefore, not be 
used anymore. 
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^> ################################################################################# 

_> # Manifestly Covariant Gauge-invariant Cosmological Perturbation Theory 
I> # 

_> # P . G .Miedema and W.A.van Leeuwen 

I> ################################################################################# 
_> # Maple 14 

I> # 

_> # Consistency check of the perturbation equations (21) . 
I> # 

_> # To check Eqs.(21) with coefficients (22), substitute (20) into Eqs.(21) 
_> # and use the background equations (5) and the perturbation equations (18) 
_> # to evaluate all time derivatives . 

I> ################################################################################# 

_> restart; 

_> # Constants : 

_> # kappa = 8*pi*G/c^4, 

_> # Delta is the nabla operator, 

'.> ################################################################################# 
_> # Coefficients b„l, b__2, b_3, (22) : 

"> bl : =kappa*e_0 (t) * (1+w (t) ) /H (t) -2*dif f (beta (t) , t) /beta (t) -H (t) * (2+6*w (t) +3*beta (t) ^2) +R_0 (t) * (1/ (3* 
H(t) )+2*H(t) * (l+3*beta(t) ^2) / (R_0 (t) +3*kappa*e_0 (t) * (l+w(t) ) ) ) ; 



hi := 



Ke_0(t) (1 +w(l)) 
Hit) 



dt 



<(.t) 



-Hit) (2+6 w(f) +3 PW^) +R_0{t) 



1 



IHW (1+3P(0^) 



3 Hit) RJKt) +3 Ke_0(fl (1 +w(l)) 



(1) 



> b2 :=-(l/2) *kappa*e_0 (t) * (l+w(t) ) * (l+3*w(t) )+H(t) "2* (l-3*w (t) +6*beta (t) *2* (2+3*w(t) ) )+6*H(t) * (diff 
(beta(t) ,t) /beta(t) ) * (w (t) +kappa*e„0 (t) * (l+w(t) ) / (R_^0 (t) +3*kappa*e_0 (t) * (l+w(t) ) ) ) -R J (t) * ( (1/2) *w 
(t)+H(t) '•2* (l+6*w(t) ) * (l+3*beta(t) ^2) / (R_0 (t) +3*kappa*e_0 (t) * (l+w(t) ) ) )-beta(t) ^2* (Delta/a (t) *2- 
(l/2)*R_0(t) ) ; 

b2:=-^ Ke_0{t) (1 +w(0) {1+3 w{t)) +mf (l -3 w(0 +6 P(fl^ (2+3 w{t))) 



(2) 



(,H(t) 



[>') 



w(t) + ■ 



Ke_0(t) (1 +w(t)) 
R_0(t)+3Ke_0(t) (1 +w(t)) 



P(/) 



2 R_0{t)+3Ke_OW a+w{t)) ) \ a{tf 



R_0(t) 



> b3:= (-18*H(t) ^2* (e„0 (t) *p_en (t) * (1+w (t) ) +2*p_n (t) *dif f (beta (t) ,t)/beta(t) / (3*H (t) ) -beta (t) ^2*p_n 
(t) +p_e (t) *p_n (t) +n_0 (t) *p_nn (t) ) / (R_0 (t) +3*kappa*e_0 (t) * (1+w (t) ) ) +p_n (t) ) * (Delta/a (t) ^2- (1/2) *R_0 
(t)) *(n_0(t)/e_0(t)); 



h3 := ■ 



1 



e 0(t) 



( i' d „ 

p_n{f) - P(0 

e 0(t) p en(t) (1 +w(t)) + ^ - P(rt^p n[t) +p e(t) p n{t) +n OU) p nn{t) 

3 p(f) Hit) 

RJKt) +3 Ke_0(f) (1 +w(f)) 



(3) 



Vp_n{t) 



a(t)' 



RJKt) nj(t) 



J 



-> ################################################################################# 
_> # Background ecjuations 

I> ################################################################################# 
'> w(t) :=p_0(t)/e_0(t) ; 



_> # Scale factor 

^> define(a,diff (a(t) ,t)=H(t) *a(t) ) ; 
'> diff (a (t) ,t) ; 



w{t) 



._ pj(t) 
ejit) 



W) a{t) 



(4) 



(5) 



> # The time derivative of Eq. (5a) combined with Eqs . (5b) and (5c) yields the time derivative of the 
Hubble function H: 
^> define (H, diff (H (t) , t) =- (1/6) *R_0 (t) - (1/2) *kappa* (e_0 (t) +p_0 (t) ) ) ; 
'> diff (H(t) ,t) ; 



(6) 



- - R_0{t) - - K (e_0(0 +p_0{t)) 
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_> # Either define the three-space curvature by Eq. (5b) : 
_> define (R_0,diff(R_0 (t) , t) =-2*H(t) *R_0 (t) ) ; 
_> # or by Eq. (6) with k=+l, 0, -1: 
^> # k:=+l; R_0 (t) :=-6*k/a(t) ^2; 
"> diff (R_0 (t) ,t) ; 



. -lH(t)RJ){t) 
_> # Energy conservation law (5c) : 

^> define (e_0, diff (e_0(t) ,t)=-3*H(t) *e_0 (t)* (l+w(t) ) ) ; 
"> diff (e_0 (t) ,t) ; 



-^H(t)e Oil) I l+PJ^ 
e_0{t) 



_> # Particle number conservation law (5e) : 
^> define (n J, diff (n_0(t) , t) =-3*H (t) *n_0 (t) ) ; 
"> diff {n_0 (t) ,t) ; 



. -3 Hit) n_Oit) 

_> # Time derivative of p„0 (t) : 

^> define (p J, diff (p J (t) , t) =p_e (t) *dif f (e„0 (t) , t) +p_n (t) *dif f (n_0 (t) ,t) ) ; 
'> diff (p_0 (t) ,t) ; 

-3 p eit) Hit) e Oit) ( 1 + EJ!^ ]-3p nit) Hit) n Oit) 
- [ ejit) J 

_> # Partial derivative p_e(t) of the pressure p: 

_> define {p_e, diff (p_e{t) ,t)=p_ee{t) *diff (e_0 (t) ,t)+p_en(t) *diff (n_0 (t) ,t) ) ; 
'> diff (p_e(t) ,t) ; 

-in eeit) Hit) e Oit) [ 1 + EJI^ 1 - 3 p en{t) Hit) n Oit) 
[ e_Oit) J 

_> # Partial derivative p_n(t) of the pressure p: 

_> define {p_n, diff (p_n(t) ,t)=p_en(t) *diff (e_0 (t) ,t)+p_nn(t) *diff (n_0 (t) ,t) ) ; 
"> diff {p_n(t) ,t) ; 



-3 p enit) Hit) e 0{t) ( 1 + EJIMI. | - 3 n nn{t) Hit) n Oit) 
. y ejit) ' 

_> # Quantity beta(t) : 

'> beta(t) :=sqrt (diff (p_0 (t) ,t)/diff (e_0 (t) ,t) ) ; 



-9 p eit) Hit) e Oit) [ 1 + EJI}^ \-9p nit) Hit) n Oit) 
- I e Oit) ) - 

P_Oit) 



m : 

Hit) e Oit) 

- [ e_Oit) 

-> ################################################################################# 
_> # First order perturbation equations 

I> ################################################################################# 
_> # First order perturbation to the pressure: 
'> p_l(t) :=p_e(t)*e_l(t)+p_n(t)*n_l(t) ; 

. pjit) ■-p_eit) ejit) +pjiit) njit) 

_> # Energy conservation law (18b) : 

'> define (e^l, diff (e_l(t) ,t)=-3*H{t) * (e„l (t) +pj (t) ) -e^O (t)*(l+w(t) ) * (theta (t) + (kappa*e_l (t) + (l/2)* 
. R_l(t))/(2*H(t)))); 
> diff (e_l (t) ,t) ; 



-3 H(t) iejit) +p_eit) eJit) +p_nit) nJit)) - eJit) f 1 + \ 

\ e_Oit) J 



, Ke lit) + - R lit) 
1 ^ 7 ^ 

6(0 + 



2 Hit) 



_> # Particle number conservation law (18d) : 

^> define (n_^l, diff (n_l(t) ,t)=-3*H(t) *n_l (t)-n_0(t) * (theta (t) + (kappa*e_l (t) + (l/2) *R_1 (t) ) / (2*H(t) ) ) ) 
"> diff {n_l (t) ,t) ; 



-3 H(t) n_l(t) -n_0(t) 



9(0 + 



, Ke ht] + - R Ht) 
1 - 2 - 



2 Hit) 
_> # Momentum conservation law (18c) : 

'> define (theta,diff(theta(t) , t)=-H(t) * (2-3*beta(t) ^2) *theta (t) -Delta/a (t) ^2*p_l (t) / (e_0 (t) * (l+w(t) ) ) 

L ) '■ 
> diff (theta(t) ,t) ; 



(16) 



-Hit) 



2 + 



, -9p e(t) Hit) e Oit) (\+ EJM. \-9p nit) H[t) n 0(f) 
1 ^- - \ e Oit) ' 



Hit) e_Oit) 1 + 



pjm_ 

ejit) 



0(0- 



A ip_eit) eJit) +p_nit) njit)) 



ait)^e_Oit) 1 + 



pjm_ 

e_Oit) 



(17) 



_> # Evolution equation for the local perturbation to the 
_> # global spatial curvature, (18a) : 

'> define (R_l, diff (R_l(t) ,t)=-2* (H(t) *R_1 (t) -kappa*e_0 (t)* (l+w(t) ) *theta (t) ) -R_0 (t)/ (3*H(t) ) * (kappa* 
. e_l(t) + (l/2)*R_l{t))); 
> diff (R_l (t) ,t) ; 

„ 1 R^Oit) Uejit) + ^ R_m 

-2Hit)R lit) +2Ke Oit) 11 + | 6(0 

- I e_Oit) ) 3 Hit) 

-> ################################################################################# 
_> # Gauge-invariant density perturbations 

I> ################################################################################# 
_> # Gauge-invariant perturbation to the energy density, (19a) : 

'> e_gi (t) : = (e_l (t) *R_0 (t) -3*e_0 (t) * (1+w (t) ) * (2*H (t) *theta (t) + (1/2) *R_1 (t) ) ) / (R_0 (t) +3*kappa*e_0 (t) * 
(l+w{t))); 



(18) 



e_giit) : 



e lit) R Oit) -3e Oit) I 1 + ] ( 2 Hit) 6(0 + - R Ht) 

- - - ' e Oit) J [ 2 - 



R Oit) +3 Ke_0(O 1 + 



pjit) 
e_Oit) 



(19) 



_> # Gauge-invariant perturbation to the particle number density, (19b) : 

'> n„gi (t) : =n_l (t) -3*n_0 (t) * (kappa*e_l (t) +2*H (t) *theta (t) + (1/2) *R_1 (t) ) / (R_0 (t) +3*kappa*e_0 (t) * (1+w 
(t))); 



n_giit) :=«__/ (0 ■ 



3 n_0(0 I Kejit) +2 Hit) 6(0 + - RJit) 



R_0it)+3 Ke_Oit) 1 + 



pjm_ 

e_Oit) 



(20) 



_> # Gauge-invariant contrast functions, (20) : 
> delta_e (t) : =e_gi (t) /e_0 (t) ; delta_n (t) : =n_gi (t) /n_0 (t) ; 



delta eit) := 



e !it)R 0(0 -3 e 0(0 fl +^=^1 f 2 //(O 6(0 + - « Aol 

~ ~ ~ I e_0(0 J \ 2 - ) 

R Oit) +3 Ke_0(0 f 1 + f^J^ 1 \ eJit) 



nJit) ■ 



e_Oit) 

3 nJit) I Ke_/(0 +2 Hit) 6(0 + ^ RJit) 



R Oit) +3 Ke 0(0 1 + 



delta_nit) := ■ 



pjm_ 

e Pit) 



n_Oit) 



(21) 



_> # Consistency check for Eq. (21b) [left-hand side minus right-hand side] : 

'> simplify (diff (delta_n (t) -delta_e (t) / (1+w (t) ) , t) -3*H (t) *n_0 (t) *p_n (t) / (e_0 (t) * (1+w (t) ) ) * (delta_n (t) 
-delta_e(t) / (l+w{t) ) ) ) ; 



(22) 



_> # Consistency check for Eq. (21a) [left-hand side minus right-hand side] : 

'> simplify (diff (delta_e (t) , t$2) +bl*dif f (delta_e (t) , t) +b2*delta_e (t) -b3* (delta_n (t) -delta_e (t) / (1+w 
(t)))) ; 



(23) 



L> 
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